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£; Abstract 

CN ' We study the density driven symmetry breaking in holographic superconductors by considering 

, the positive mass squared case. We show that even for the positive m^, a scalar condensation still 

*0 ' forms, provided the chemical potential is high enough. As increases, the phase space folds due to 

'NT ' 

' the non-linearity of the equations of motion, and two nearby points in the phase space can represent 

o : 

. symmetry breaking and preserving configurations respectively. The phase space defined by the set 

of initial conditions of field variables at the horizon undergoes a non-linear radial evolution to result 
^ ' in the phase space folding, a characteristic phenomenon in a non-linear system. We then calculate 

the specific heat, which characterizes superconductors and has been measured in experiments. We 
observe a discontinuity in the specific heat at the transition point and compare our results with 
experimentally observed numbers. The electrical conductivity for various is also calculated. 
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I. INTRODUCTION 



It has been suggested that some of phenomena of strongly interacting systems can be 
described by the AdS/CFT correspondence While the relevance of the gauge/string 
duality to the QCD has been much discussed in the context of heavy ion collisions 2, ^, 4| 
and hadron physics, its relevance to the electron system is rather new and still an open 
issue. However, testing the applicability with a model in the bulk and the phenomena at 
the boundary is a very exciting possibility. In this context, it has been attempted to use 
the AdS/CFT correspondence to describe certain condensed matter systems such as the 
quantum Hall effect 5||, Nernst effect , 0], superconductivity [o], [lol, [ll| , and fractional 
quantum Hall effect 121 1. 

ri 

In Ref. [9|, a model of a strongly coupled system which develops superconductivity was 
suggested based on the holography, which is an Abelian Higgs model in 3+1 dimensional 
anti-de Sitter(AdS) spacetime. In this model a complex scalar field, as a dual of the "Cooper 
pair" operator at the boundary, is introduced whose vacuum expectation value is the order 
parameter of the superconductivity. The symmetry breaking was induced by considering 
the negative mass squared of the scalar field. In AdS spacetime, the negative mass squared 
is allowed due to the Breitenlohner-Freedman(BF) bound. In this case, the dual operator 
coupled to the scalar field at the boundary has a weight A = 1 or A = 2. Recall that for 
a scalar field in AdS^+i spacetime, (mL)^ = A(A — d). The cases with A = 3/2 and 3 
were considered in Ref. 11[. The hydrodynamic nature was studied in a theory with large 
conformal operators, 1/2 < A < 5, in Ref. [l^. 

In 2+1 dimensions, a fermion bilinear (Cooper pair) has naive dimension 2, therefore its 
dual scalar field should have a tachyonic mass, leading to the instability for the symmetry 
breaking. However, the anomalous dimension is usually non-zero and can be large, so we 
are forced to ask what happens if the mass squared is positive. Another motivation is the 
neutron superfluidity in the neutron star. Since the neutron consists of three quarks, its 
creation operator is a composite operator with high conformal dimensions and so is the 
Cooper pair of the neutrons, naively A = 9 for the pair. Again, one is led to the question 
of the positive mass squared case. 

When the scalar field has a positive mass squared, one may expect that the mechanism 
of the symmetry breaking is also removed. However, in the presence of the matter, there 
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is a density driven instability leading to the symmetry breaking as the effective mass of 
the complex boson in flat spacetime is given by m^j^jr = m? — /x^. Here /x is the chemical 
potential. 

In warped spacetime, we can still expect the same mechanism at work. This is already 
pointed out in Ref. 19|] from the fact that the fluctuation of the scalar field has an effective 
mass 

m%f = m^-W'\Al (1) 

where At is the electrostatic potential whose value at the boundary is the chemical potential. 
Notice that the regularity requires = at the horizon. As a consequence, the second 
term of Eq. ([T]) goes to zero at the horizon in spite of the divergence of On the other 
hand, near the boundary, is zero hence it eliminates the effect of the chemical potential. 
Therefore, the density effect should be less effective than the flat space case, and it is not 
clear whether the density driven symmetry breaking mechanism is still working. The present 
work is to address this question in detail. 

In fact, our results show that for a small positive mass s quar ed, the result is not much 



different from the negative case studied in Refs. , ll5|, llo, [ITI, ll9|] . However, as the mass 
squared increases, the phase space folds so that two very close boundary conditions can give 
completely opposite states: one symmetry broken state, and another symmetric state. 

The specific heat is one of the experimental observables for superconductors, and it has 
two interesting limits. At a temperature far below the transition temperature Tc, at which 
the superconductivity disappears, the specific heat is exponentially small due to the pairing 
gap, while it exhibits a discontinuity very near T^. We calculate the specific heat with various 
m? to see if large conformal dimensions change the properties, focusing on the normalized 
magnitude of the discontinuity. Then we compare our results with experimentally observed 
numbers of some real materials. Finally we calculate the electrical conductivity for m? > 
and < 0. 
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II. ABELIAN HIGGS MODEL AND HOLOGRAPHIC SUPERCONDUCTORS 



We use the Abelian Higgs model in AdS^+i Schwarzschild black hole background to study 
the holographic superconductors. The action of the model is given by 

S=- j d'^^^x^ {^F'^^Fab + |<9^ - zA^I' + m^l^l'^ , (2) 

where the bulk Maxwell field strength is Fab = da^b ~ db^ai and is a charged complex 
scalar field. We work with the following (d+l)-dimensional background metric: 

ds' = ^^{-f{z)de + ^^+dx'^^ f{z) = l-{£j, (3) 



where the horizon is at 2; = 2:0. The Hawking temperature of the black hole is given as 

d 



(4) 



A-KZq 

Note that we will work in the probe limit: the Maxwell field and the scalar field do not give 
any back-reaction to the background metric. Here we take the form of the bulk scalar field 
as = -^ip^^'^ and consider the gauge transformation Aa Aa + da4>- Then the equations 
of motion for static and translational invariant fields of interest read 

//' d-l\ /A? ^ 



/ z • \P z^f 

A': + i3-dA-^A = 0. (5) 
z z^f 

To solve the equations of motion, we specify appropriate boundary conditions as follows: at 
the horizon, 

A{zo) = 0, i:{zo) = -^i:'{zo), (6) 
and at the boundary, i.e. z ^ 0, 

A{z)=ix-pz''-^ + --- . (7) 

We will fix the chemical potential fi and thus work in the grand canonical ensemble. For 
the boundary value of the scalar field ip{z), we will choose = 0. When both and 
are normalizable modes we need to make a choice of either \E'i = or \1'2 = ^j. As we 



increase m^, A increases and becomes non-normalizable mode, which corresponds to a 
source term. We choose = to impose the source free condition which is necessary to 
study the spontaneous symmetry breaking rather than induced breaking. 

Notice that the equations of motion Eq. dS]) are invariant under the re-scahng: 



- - 4 - 

Z — ZqZ, /if — , 



with tp{z) fixed. At the boundary, we then have 



(8) 



At{z) fl- pZ^-^ 



(9) 



where 



/i = /iZo, 



p = p4 \ 



^1 = ^1^-^, and ^5 



(10) 



There have been many studies on t 



2+1 and 3+1 dimensions 
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le properties of the holographic superconductors in 
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19| . We note that in the hteratures, im? has been 



assumed to be negative to guarantee the symmetry breaking and the scalar field condensate. 



In Ref. 



Jjj, the massless case was also considered. 



III. THE PHASE SPACE FOLDING 



In this section, we extend the previous studies on the holographic superconductors by 
considering the bulk scalar field with the positive mass, m? > 0. As in the m? < cases, we 
need to solve the coupled equations of motion in Eq. ^ with the boundary conditions, (jS]) 
and (IT!) with the choice of = and fixed /i. 

First we define the 'phase space' as a set of possible values of pair (?/;', A[) at the horizon 
Zq, which evolves radially to the boundary 2 = 0. To understand this properly, it is useful to 
consider some details of the numerical integration. If we choose to integrate starting from 
the horizon to the boundary, we need the values of the fields and their first derivatives at 
the horizon. However, half of the actual boundary conditions are imposed at the boundary, 
z = Therefore we artificially assign arbitrary values of the first derivative of the fields 
{iI)\zq)^ A'^{zq)) at the horizon. Some of them will give the correct boundary conditions 
given in Eq.d?]) with \E'i = at the boundary. Such a set of 'initial conditions' giving 
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FIG. 1: Phase space as the contours of {zq), A'^{zQ)) satisfying ^'i = 0. (a) m? = —2, (b) 
rm? = 0. Each contour hue satisfies the boundary condition ^'i = 0. 

correct boundary conditions draws contours in the plane of [ip^zo)^ ^^{zq)). We define such 
contours as our phase space. We use the word "phase space" to have the same meaning as 
the configuration space. 

For m? < 0, the phase space is drawn in Fig. [1] (a). In Figs. [T]and[2], the contours are for 
different nodes. As A'^{zq) approaches zero, some contour lines misleadingly look connected 
to its neighbor but none of them actually meet each other. 

As the value of m? increases, the phase space begins to wrap, as seen in Fig. [21 In order 
to see what is happening, we choose a phase space for m? = 1.36 and take several points 
there as seen in Fig. [31 We number the phase space points such that fi/T increases as the 
number increases. Then we calculate the scalar field condensate and the specific heat and 
draw the shape of the scalar field configuration. See Figs. [3l [H and [51 respectively. 

We see that the wrapped part which contains the points 6-10 is unstable since the specific 
heat is negative there. To support this claim further, we also study the shape of the scalar 
field configuration and draw the result in Fig. [51 For a point in the folded region of the 
phase space, the corresponding wave function develops a needle like shape which costs too 
much energy, making the configuration unrealizable. 

On the other hand, the configurations in the unwrapped part (containing points labeled 
by 1,2,3,4,5) are regular and realizable. As we increase further, however, the phase space 
progresses to fold itself until it gives an interesting behavior: two infinitesimally close points 
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FIG. 2: Contours of {ip' (zq), A[{zq)) for various positive m^. Each contour line satisfies the 
boundary condition = 0. (a) = 0.31, (b) = 1.36, (c) = 3.04, (d) = 4, 

in the phase space wih give completely different characters of the states in the boundary, 
one ordered and the other disordered. At the onset of such an effect, the superconducting 
and the normal states might coexist in an arbitrarily mixed fashion. It would be interesting 
if such a mixing is in chaotic fashion. 

Discarding the disordered part, the scalar field condensates for various are shown in 
Fig. [61 In Fig. [71 we plot the critical temperature as a function of the mass squared. T^. 
decreases as increases, as it should be. 
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FIG. 3: The wrapped phase space for m? = 1.36 and the condensation. The numbers, 1-10, were 
to indicate the correspondence between the configuration and the condensation state. Since Tc 
is fixed by the chemical potential, the marking numbers can also be considered as a measure of 
(increasing) chemical potential for fixed temperature. 




0.0 0.2 0.4 0.6 0.8 1.0 

T_ 
Tc 

FIG. 4: Specific heat for each point on the phase space in Fig. [3l The numerical values of the 
specific heat for the points 7, 8, 9 and 10 are -1389, -6106, -5581 and -4731 respectively. 

IV. SPECIFIC HEAT AND CONDUCTIVITY 

In the previous section, we studied the density driven symmetry breaking in particular 
with large conformal dimensions. We now consider two observables that exhibit characteris- 
tics of superconductors with negative and positive values of to see how those observables 
are affected by varying or conformal dimensions. We first consider the specific heat of 
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0.002 0.004 0.006 0.008 0.010 



FIG. 5: Shape of for a few points of phase space in Fig. [3l The lower right figure is the closer 
look near z = of the needle shape in the lower left figure. 



<>Pj>l/A 




FIG. 6: The condensation ^'2 in d=2+l (left) and d=3+l (right). Left: the thick solid line, the 
thin solid line, and the dashed line indicate = —2,m? = 0.31, and m? = 1.36, respectively. 
Right: the thick solid line and the thin solid line correspond to the cases = —3 and = 1.29. 



holographic superconductors. It is known that the specific heat of superconductors shows 
an interesting behavior at low and high temperatures: at a temperature far below the tran- 
sition temperature Tc, the specific heat is exponentially small due to the pairing gap, while 
it exhibits a discontinuity very near Tc. To calculate the specific heat, let us compute the 
free energy of the boundary field theory. 
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FIG. 7: The critical temperature Tc as a function of in d=2+l (left) and in d=3+l (right). 



The free energy is obtained by the on-shell action of Eq. ([2]) 

d = T Sqs 

where 



-99 



(11) 



;i2) 



2 

Here e is introduced to cut off the divergence at the boundary that will be regularized by a 
boundary counter-term. With the regularized on-shell action, we obtain the free energy of 
the superconducting phase 



^s(/i) 

V 



[2-d)fip + 



20 



Af^^ dz 



(13) 



where V is the volume factor. 

From the free energy, we obtain the specific heat as 



(14) 



V (9T2 ■ 

We note that in experiment, one generally measures the specific heat at constant pressure, 
Cp. The difference between Cp and c^, the specific heat at fixed volume, is ~ {dV/dTy. 
Therefore we may ignore the difference in condensed matter systems since the volume of the 
system will hardly change with respect to the temperature. In Fig. [HI we plot our results of 
the specific heat for the holographic superconductors with a few masses. 

For usual superconductors, it is known that the specific heat vanishes exponentially as 
the temperature approaches zero; 



(15) 
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FIG. 8: Specific heat as a function of temperature in d=2+l(left) and in d=3+l(riglit). Left : 
Tlie tliick solid line, the thin solid line, the dashed line are for m? = —2,m? = 0.31, = 1.36 



where Aq is the gap at zero temperature. In our evaluation of the specific heat, we do not 
observe the exponential suppression near T = 0. Due to some numerical instability near 
zero temperature, it is not clear whether this is a genuine feature of the holographic super- 
conductors or numerical error. We notice that a similar observation has made previously in 



Now, we consider the ratio of the specific heat of the normal phase (c„) to the supercon- 
ducting phase at the critical temperature. Here, normal phase is defined as the one with no 
condensation, i.e., ip{z) = 0. The ratio is useful since there are experimentally measured 
values of it, and also near the critical temperature, we can trust our numerical calculations. 
The ratio is defined by 

R = "-^^ , (16) 

Cn T=Tc 

which measures the discontinuity of the specific heat at the critical temperature. The cal- 
culated ratios are listed in Table 1. 

Table 1 : The values of R for various for the holographic superconductors in 



respectively. Right : The thin solid line is for 



3, and the thick solid line is for = 1.29. 




d=3+l. 
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R 


-3 


3.01 


1.29 


2.27 


1.76 


2.24 


3.29 


2.17 


5 


2.11 



To see if our results have any relevance to a real condensed matter system, we may 
compare our numbers with the ones from experiments. In Table 2, we list the values of R 
for some materials 22 1. 



Table 2 : The values of R for various materials. 



Material 


R 


Material 


R 


BCS theory 


1.43 


Al 


1.45 


Cd 


1.36 


Pb 


2.71 


CeRu2Si2 


3.5 


BaPbi_,Bi,03 


8 


(Lao.925Cao.075)2Cu04 


5.7 


YBasCugOr 


3.6 



Finally, we consider the electrical conductivity, which shows the characteristics of super- 
conductors. The conductivity can be obtained by taking into account the retarded Green's 
function, <j{uj) = ■^G^{uj,k = 0). Now we introduce a gauge field in the bulk as a 
fluctuation which couples to an external electric current at the boundary. The equation of 
motion for Ax is 



A'' + 



f 3-d 



a: 



Ax = 0. 



To solve the equation, we impose the ingoing boundary condition, A^ 
horizon. At the boundary, we have 

Ax = A^'^+A^'h + --- . 

Then the retarded Green's function reads 

and thereby 



A^ 



UOA: 



(0) 



|fc=0 



(17) 

(18) 
(19) 

(20) 
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FIG. 9: The conductivity at T ~ 0.15 Tc for = —2 (left), and m? = 1.36 (right). The sohd hne 
is for the real part, while the dashed line represents the imaginary part. 

In Fig. [9l we present conductivities for = —2, and = 1.36 in 2+1 dimensions. We do 
not find any qualitative difference between the positive and negative cases. Our results 
are qualitatively similar to the one obtained previously with < UJ, ll7|]. The peak of 
the conductivity becomes a little bit broad with increasing m^, but the dependence is weak. 



V. SUMMARY AND DISCUSSION 



We studied the holographic superconductors by the density driven instability by consid- 
ering bulk scalar fields with > 0. We found that the holographic superconductors still 
exist with positive m^. We observed that as the bulk mass of the scalar field increases, the 
system exhibits the phase space folding: two very close boundary conditions can give com- 
pletely opposite states; one symmetry broken state, and the other symmetric state. Such 
sharp dependence on boundary conditions is reminiscent of "the butterfly effect" , a road to 
chaos. As a consequence, for large conformal operators, the initial conditions prepared at 
the horizon should be highly fine tuned to obtain the density driven symmetry breaking. 
However, we remark here that while the usual butterfly effect is in the time direction, our 
case is for a radial evolution in the bulk of the dual gravity. It would be interesting to study 



it further in relation to the chaos 



2lj. 



We also studied other quantitative properties of the holographic superconductors, i.e. the 
specific heat and the conductivity. We found that the value of the specific heat increases 
rapidly as temperature approaches the critical temperature, which agrees with a genuine 
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feature of superconductivity in real materials. The ratio of the specific heat of the normal 
state to the superconducting state at the critical temperature are compared with experi- 
mentally observed ones to see if holographic superconductors share some common feature 
with real materials. The obtained behavior of the conductivity was rather universal against 
varying m^. 

Finally, we remark that to establish the instability of the folded region in the phase space 
more reliably, one needs to study the back reaction of the scalar field to the geometry. The 
instability for the large chemical potential sectors might be due to the fixed metric. In this 
work, the charge is allowed to fluctuate. If the scalar field were absent, we know that the 
charge is not supposed to be larger than the critical value when considering a fixed mass. 



23| 



Here, the coupling of the gauge field and the scalar field makes things fuzzy. The work 
claims that for any values of the scalar charge larger than a bound value, the condensate 
is always possible. It would be interesting to study how our work is related to it, and the 
work in this direction is in preparation. 
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